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RESUMO 


A formulação matricial da análise estrutural veio 
remediar alguns dos graves inconvenientes que caracte- 
risavam a Teoria das Estruturas. Só porém a técnica dos 
elementos finitos viria dar à análise estrutural a genera- 
lidade de que carecia. 

No presente trabalho, concebido para alunos, apresen- 
tam-se em primeiro lugar elementos de análise matricial 
e, seguidamente, noções elementares da técnica dos ele- 
mentos finitos. 


1— INTRODUÇÃO 


A Teoria das Estruturas antiga era caracteri- 
zada por gravíssimos defeitos. Destacam-se entre 
eles: 


— a insuficiente generalidade (só considerava 
estruturas reticuladas) 

— a deficiente sistematização e falta de unidade 

— o isolamento relativamente aos métodos ge- 
rais da Análise Numérica 

— a insuficiência dos instrumentos matemáti- 
cos utilizados 


A formulação matricial da teoria das estrutu- 
ras reticuladas veio remediar todos estes incon- 
venientes excepto o primeiro. Só a técnica dos 
elementos finitos daria mais tarde à Teoria das 
Estruturas a generalidade de que carecia. 


SYNOPSIS 


The matrix formulation of structural analysis could 
obviate to some defects of the Theory of Structures, 
Only the finite element technique would give however 
to structural analysis the generality which was missing. 

In the present paper intended for students, basic 
notions of matrix analysis are first introduced, followed 
by a short description of the finite element method, 


Não pode surpreender que a Álgebra Matricial 
seja um instrumento adequado ao cálculo de es- 
truturas. A arrumação do cálculo mediante o 
preenchimento de tabelas que se obtêm suces- 
sivamente uma das outras é com efeito prática 
corrente nos gabinetes de estudo, e a Álgebra 
Matricial não é no fundo mais que uma álgebra 
de tabelas. 


2— ANÁLISE MATRICIAL DE ESTRUTURAS 
RETICULADAS 


2.1 — Matriz de rigidez 


Considere-se uma barra no espaço solicitada 
por forças e momentos nas suas extremidades. 
O estado de tensão e deformação no interior 
da barra fica determinado se forem dadas essas 


(*) Notas de lições proferidas pelo autor na Universidade Federal do Rio de Janeiro em Agosto e Setembro de 1969. 
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forças e momentos, ou os deslocamentos e rota- 
ções correspondentes. 

Designar-se-ão as forças e momentos por forças 
generalizadas, e os deslocamentos e rotações por 
deslocamentos generalizados. 

Sejam Fº, Mº e Fº Mº os vectores das forças 
generalizadas actuantes respectivamente nas extre- 
midades esquerda e direita. Sejam u",t“ e us e 
os correspondentes vectores dos deslocamentos e 
rotações (fig. 1). 


pó ,FS 
nt Fº b as MS 


Fig. 1 


é e. e. e se ne 
Sejam u,, U,, U; as componentes deu TEM os Os 
as de 1º, etc., e considerem-se os seguintes 
vectores; 


q"º = ut | ci us 
Us us 
us us 
os 0% | 
o, og | 
9€ q 
(3) (3) 
qi = | E! qi | Fa | 
d 

| FS ES 

Mi M' 

M$ M$ 

e d 
[8] (Ms) 

e ainda 

be ft nbe 

q" = : | Q" — : 
q"! ] qa 


Supondo que os deslocamentos generalizados 
das extremidades da barra sofrem uma variação 
infinitesimal, a correspondente variação da ener- 
gia da deformação é igual ao trabalho das forças 


exteriores pela variação dos deslocamentos. 
Tem-se assim 


Gu" = 0" aq (1) 


Ora, pelo princípio da sobreposição, as forças 
são funções lineares dos deslocamentos, isto é, 


ce tkig +... 


+kygi+. x ia 


Q,=kyq + En go + - 


Q=45 a!+ "o 


(2) 
em que se denotam por Q; ADO ua « « GEES 
b bob 
mentos do vector À, e por g,, q,, ... os ele- 
mentos do vector q”. 
Fazendo 
tem-se 
Qº=kKº q? (4) 


Kº diz-se matriz de rigidez. 

Suponham-se nulos todos os elementos de q” 
excepto o de ordem j que se toma igual a um 
(estado de deformação unitária de ordem j). Vem 


— RR (4º) 
Como ki; ,..., kiy ,... são os elementos da coluna 
de ordem j de K', segue-se que os elementos da 
coluna de ordem j de Kº são as forças que é neces- 
sário aplicar na barra para conseguir o estado de 
deformação unitária de ordem j. O elemento ka | de 
K” é pois a força de ordem i necessária para AA 
seguir o estado de deformação unitária de ordem j. 

Os teoremas da reciprocidade permitem concluir 
que a força de ordem i necessária para conseguir 
o estado de ordem j é igual à força de ordem j 
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necessária para conseguir o estado de ordem i. 
Por isso, 


b b 
kj = k, (5) 
Isto signica que K? é uma matriz simétrica, 


isto é, que 
pT 


Kº =K (57) 
Introduza-se (4) em (1). Vem 
UP = (x q") sq" a” q?” KT q? = 
= q Kºôq (6) 
Integrando, obtém-se 
U" =p q” K q” (7) 


A matriz de rigidez é assim a forma quadrá- 
tica da energia de deformação. 


2.2 — Caso plano 


Considere-se uma barra plana solicitada e 
deformada no próprio plano. 

Tem-se em cada extremidade duas componentes 
de força e um momento, e duas translações e 
uma rotação. O vector dos deslocamentos é 


Us E, 
g e Mº | 
ul | sf 
us a] 
94 | E] 
E y 4 


Suponha-se a barra de eixo rectilíneo e secção 
constante (área Q& e momento de inércia 1). 
A matriz de rigidez, é, no caso de se desprezar 
a deformação por esforço transverso, e de os 
sentidos positivos das componentes serem os 
indicados na figura 2, 


* 
1 


A primeira coluna da matriz de rigidez con- 
tém as forças necessárias para que a extre- 
midade esquerda sofra uma translação unitária, 
sendo nulos todos os restantes deslocamentos 


(fig. 3). 


E 2 


L ain 


. ER 
L 


Fig. 3 


A quinta, por exemplo, refere-se a uma trans- 
lação, unitária e perpendicular à barra, da res- 
pectiva extremidade direita (fig. 4). A sexta, à 
rotação unitária da mesma extremidade (fig. 5), 
e assim sucessivamente. 


Fig. 4 


Fig. 2 
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No caso de as componentes não se referirem 
a eixos que não sejam paralelos e perpendiculares 
à barra, a situação complica-se. A matriz referida 
a esses eixos pode no entanto ser determinada, 
por operações matriciais simples, a partir da 
matriz referida ao referencial próprio da barra, 
isto é, ao eixo da barra e a eixos perpendiculares. 


2.3 — Mudança de coordenadas 


Fig. 6 


Considere-se uma rotação O dos eixos coorde- 
nados. 


Sendo v um vector arbitrário, tem-se 


| v,=v, cos!—v, sen) 


a (9) 
Vy==V, sen J + vo cos 4 
ou, sob forma matricial, 
) — 9 
— | cos sen |v=sv (10) 
“sen? cos 8 
em que 
| Vo | Va 
É fácil reconhecer que 
9TS$S=S8F= (11) 


Por isso, multiplicando (10) à esquerda por e, 
vem 
S yv=S8"Syv=v' (12) 
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Consegue-se, pois, com toda a facilidade, expri- 
mir ví em termos de v. 

No espaço tridimensional obter-se-iam expres- 
sões semelhantes às que foram obtidas para 
o caso bidimensional. 


bb 
Atendendo a que u, ,u,,-- 
de vectores, tem-se 


são componentes 


$º0 0 0 
be 0 1 0 0 be b be 
— 4 R +, 
À o 0 so | à 
0001. 
que =P" q” (13) 
bd pb gba gra E Rb" q" (14) 
(Pe Rº 0 be E be aaa RP! “tha (15) 
qi — Rº qa qa — RT qa (16) 
Multiplicando (4) à esquerda por 
“pbT bt 
h e had a?) 
o Rº 
vem 
ou 
QT T q" (19) 


A nova matriz de rigidez, isto é, a matriz tal 
que 


q at K? q” (20) 
é pois 


KP = T' KPT» (21) 


2.4 — Análise de uma estrutura 


O problema que se põe agora é o da análise 
de uma estrutura reticulada, isto é, de um sistema 
de barras. 

A deformação desse sistema pode definir-se 
dando os deslocamentos (generalizados) dos nós 
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ou as forças (generalizadas) actuantes nesses nós 


(fig. 7). 


q- 19; 


md 


dz 


qa 


a 
Q, 
Q,; 
Q, | 4 


q; — vector dos deslocamentos do nó i 
Q, — vector das forças actuantes no nó i 


Fig. 7 


Em virtude do princípio de sobreposição, Q 
relaciona-se com ( por 
Q=Kq (22) 
Em virtude de raciocínio semelhante ao que 
foi feito para a barra isolada, pode concluir-se 
que K é a matriz simétrica da forma quadrática 
da energia de deformação da estrutura, expressa 
em termos dos deslocamentos nodais. 
Tem-se assim 
U=4p q" Kg (22º) 
A análise da estrutura pode ser conseguida a 
partir de três grupos de equações: as equações 
de compatibilidade, as de equilíbrio e as relações 
forças-deslocamentos das barras (4). 
Supor-se-á, na análise que se segue, que as 
forças exteriores estão aplicadas nos nós. Ver- 
-Sse-á depois como tratar casos em que as forças 
exteriores se apliquem nas barras. 


2.5 — Equações de compatibilidade 


As equações de compatibilidade afirmam que 
as extremidades coincidentes de todas as barras 
que convergem num determinado nó têm deslo- 
camentos iguais. 

Considere-se o exemplo da fig. 7. Tem-se 


l | | | 


= 
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À NM —> pl — o (24) 
nd ss | 
| q. lo 
le — 
j de (25) 
id 1 
q = i Qu 
(23) pode tomar a forma 
1000 | 
= | (26) 
q | 0100 | , 


ou, atendendo a (13) e (14), 


Io 00 
=| 5 Rio o ja=a!a (27) 


Do mesmo modo, 


T 
qe t+ Qt | q=A"q (28) 


Para uma estrutura arbitrária, as equações têm 
a mesma forma, isto é, tem-se em geral, para 
cada barra 


qº=A"q (30) 


2.6 — Equações de equilibrio 


As equações de equilíbrio exprimem o equili- 
brio de cada nó. 
Considere-se a barra da fig. 8, 


Fig. 8 
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e tome-se um troço muito pequeno junto, por 
exemplo, da extremidade esquerda (fig. 9). Seja 
e“ a extremidade do troço que não é extremidade 
da barra 


Dado que o troço ee' é muito pequeno, é 
evidente que as forças generalizadas que é neces- 
sário aplicar na extremidade e” para substituir a 
acção da porção restante da barra, têm de ser 
iguais e de sinal contrário às forças actuantes na 
extremidade e. 

Isole-se agora um nó, por exemplo, o nó 2 da 
estrutura representada na fig. 7. O raciocínio 
precedente justifica que se considerem para a 
análise do equilíbrio, as forças indicadas na fig. 10. 


-ple 


ad Mie 


pia / Si 


Fig. 10 


Tem-se pois a seguinte condição de equilíbrio : 


E A ia E 6=0 (31) 


ou 


nt" qi RI qle= q, (32) 
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As equações de equilíbrio relativas à estrutura 
tomada para exemplo são, na sua totalidade, as 
seguintes: 


HT 


Rr q! pl que q 


(33) 


(34) 


sendo as matrizes À” as que intervêm nas equa- 
ções de compatibilidade. 

2.7 — Equações de equilíbrio em geral 

A forma (34) das equações de equilíbrio não 
é acidental. Pode ser facilmente estabelecida para 
o caso geral considerando que uma variação da 
energia de deformação da estrutura é a soma das 


energias de deformação de todas as barras. 
Tem-se assim 


a " E” 
U=0" a=2su"=2o” 5q (35) 
b b 
Introduzindo (30), vem pois 
0 =(Do" a): 
du ) 90 (36) 
Ma | 


Como 24 é arbitrário, a igualdade (36) implica 
(34). 


2.8 — Resumo das equações. Método dos des- 
locamentos, 


As equações que governam o sistema são pois: 


equações de equilíbrio: Q = b> dd 0º (37) 
b 
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b 


equações de compatibilidade: q =A” q 


(38) 


relações forças-deslocamentos : q = Kº q (39) 


Para resolver estas equações pode-se começar 
por eliminar as forças, formulando o problema 
em tempos dos deslocamentos. Tem-se então o 
método dos deslocamentos. 

Eliminando os deslocamentos, tem-sé o pro- 
blema formulado em termos de forças : é o método 
de forças. 

Não se descreverá o método das forças. 

Quanto ao método dos deslocamentos, é fácil 
de formular. As forças eliminam-se introduzindo 
(39) e (38) em (37). Vem 

a= Da KA (40) 
b 
equação que serve para determinar os desloca- 
mentos. 


Comparando (40) com (22) conclui-se 
K = pa CP er (41) 
b 

Conseguiu-se pois exprimir a matriz de rigidez 
da estrutura em termos das dos elementos. Tal é 
um dos notáveis resultados da análise matricial 

das estruturas. 
Observa-se que a fórmula (41) pode tomar a 


mesma forma que a fórmula (21). Com efeito, 
fazendo 


(42) 


(a) 


vem 
q=44 (43) 
Fazendo 
n K'! cw o. 1 
K = o KI 
a (44) 
0 Kº 
a fórmula (41) transforma-se em 
Kk=4"TK A (45) 


a qual pode interpretar-se como a fórmula de 
transformação da matriz de rigidez para a 
mudança de coordenadas (deslocamentos) (43). 


2.9 — Forças actuantes nas barras 


Supôs-se até aqui que as forças exteriores se 
aplicavam nos nós da estrutura, 

O caso em que as forças estejam aplicadas nas 
barras pode converter-se facilmente num outro 
em que as forças estejam aplicadas nos nós. 

Atenda-se à fig. 11. 

Pode considerar-se o caso (a) como sobrepo- 
sição do caso (b) e do caso (c). No caso (b), os 
nós da estrutura não se movem nem rodam, pelo 
que há que aplicar nas extremidades da barra 
carregada forças e momentos de fixação iguais 
às forças e momentos de encastramento. No 
caso (b), os esforços são pois nulos em todas as 
barras excepto na barra carregada, que funciona 
como isolada e perfeitamente encastrada. No 
caso (c) forças e momentos iguais aos de fixação 
como sinal trocado são aplicados nos nós. É evi- 
dente que os casos (b) e (c) sobrepostos dão 
o caso (a). 


nós imóveis 


Fig. 11 


TECNICA Nº 395 


225 


3. GENERALIZAÇÃO PARA ESTRUTURAS 
NÃO RETICULADAS 


3.1 — À necessidade da discretização 


A teoria das peças lineares é uma teoria 
discreta, isto é, o campo elástico no interior de 
uma barra pode exprimir-se exactamente em termos 
de um número finito de parâmetros: os deslo- 
camentos generalizados ou as forças generali- 
zadas nas extremidades da barra. 

Já a teoria das cascas, por exemplo, não é uma 
teoria discreta, visto que não há possibilidade de 
definir em termos de um número finito de pará- 
metros a distribuição das grandezas que nela 
intervêm. Na realidade, a teoria das peças linea- 
res é, entre os modelos analíticos da Teoria das 
Estruturas, o único discreto. E isso explica que 
as estruturas reticuladas, que são constituídas por 
elementos lineares, fossem as primeiras a ser 
abordadas pela Teoria das Estruturas. 

Para analisar estruturas não reticuladas, há que 
recorrer à integração formal de equações diferen- 
ciais, ou introduzir hipóteses simplificativas que 
permitam transformá-las em sistemas discretos. 

O método dos elementos finitos não é mais 
que uma técnica geral de discretização. 


3.2 — Técnicas de discretização 


As técnicas de discretização utilizadas na 
análise estrutural são fundamentalmente de dois 
tipos: 


a) os métodos de diferenças finitas 
b) os métodos de sobreposição 


Os segundos baseiam-se na expressão aproxi- 
mada do campo como combinação linear de um 
número (finito) de campos elementares conhecidos. 
As incógnitas a determinar são os coeficientes 
dessa combinação linear. 

O método dos elementos finitos é um método 
de sobreposição que considera o domínio corres- 
pondente à estrutura subdividido em certo número 
de subdomínios. Os campos elementares têm 
expressão analítica diferente em cada um dos 
subdomínios, em vez de uma única expressão 
válida em todo o domínio. Diz-se que os campos 
têm definição seccional. 
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Para se compreender a vantagem da definição 
seccional, considere-se o caso de uma função de 
uma única variável y == f(x) (fig. 12). 


a a, a, aa b x 


Fig. 12 


Será em geral difícil aproximar esta função no 
intervalo (a, b) por outra função a que corres- 
ponda uma expressão analítica válida em todo o 
intervalo. 

No entanto, considerando o intervalo subdivi- 
dido em subintervalos, a função poderá ser razoa- 
velmente aproximada, por exemplo, por uma 
poligonal, isto é, por uma função cuja expressão 
analítica mude de subintervalo para subinter- 
valo (fig. 12). 


3.3 — Descrição do método dos elementos 
finitos 


O método dos elementos finitos consta dos 
seguintes passos: 


a) subdivisão da estrutura em subdomínios 
(elementos finitos) 


b) discretização do campo dentro de cada 
elemento 


c) compatibilização dos elementos entre si. 


Os programas de cálculo automático que per- 
mitem a análise da estrutura pelo método dos 
elementos finitos fixam uma forma genérica para 
os elementos (triangular, rectangular, etc.), e o 
modo como os elementos são compatibilizados. 

O utilizador que disponha de um certo pro- 
grama não tem pois mais que subdividir o domi- 
nio em subdomínios, cuja forma seja consentida 
pelo programa, e definir as solicitações exteriores. 
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3.4 — Exemplificação para a elasticidade plana 


3.4.1 — Equações da elasticidade plana 


O método pode ser aplicado na análise de 
qualquer tipo de estrutura. Será porém ilustrado 
com uma aplicação à elasticidade plana. 

As equações da elasticidade plana são as 
seguintes: 


— equações de equilíbrio : 


do do 
| T Bj X=0 
| 0 X| OX, 
(46) 
G dg 
[o q Mg 
| dx, “a E 
— relações deformações — deslocamentos 
du, 
H € 
0X, 
du 
| “gg = i (47) 
7 Xo 
| , du; du, 
“2 
| ) Xo d Xj 


— relações tensões — deformações (tensões pla- 
nas): 


E eua) 
= — E VE 
O e É 29 
2 Ge (48) 
É (E VE 
É qm) 99784) 
Tp = 2 Ges, 


Estas equações podem escrever-se sob uma 
forma simbólica muito conveniente desde que 
se considerem os vectores 


x | 74 | En u, 
Á-= | g= 1» 6= — | | 
X, | é | ge tg 
219 12 
(49-52) 
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Atendendo a que 4,, = 0,,, tem-se então, efec- 
tivamente, em vez de (46), (47) e (48), 


EE=X (53) 
e=Du (54) 
s=He (55) 


em que E e D são matrizes de operadores dife- 
renciais, 


Ea d 
0X, O Xo 
O E e 
7 Xy d x, 
[ J o ] 
dx, | 
Ô 
Dp=7, 0 Fra , (57) 
PL 
| 0X 0X| 
e H é a matriz simétrica 
| 1 v O | 
u 1 [à] 
À ss 2 CS] (58) 
E.» | é O Es 
E 


Na fronteira, as equações de equilíbrio são, 
como se sabe, 


Cn; +O%YWn=7, 


(59) 
90, + Sony = 0, 

ou, sob forma simbólica 

D==N 8 (60) 
sendo 

G 

p= fai) (61) 

e 


| n O n 
N=/ à 62 
Ê No | ea 


O interesse da notação introduzida reside sobre- 
tudo na sua generalidade. É com efeito possível 
escrever sob esta forma as equações de todos 
os modelos analíticos da Teoria de Estruturas. 
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3.4.2 — Expressão da matriz de rigidez 


O tipo de um elemento define-se pela sua forma 
genérica e pela família dos campos admissíveis no 
interior do elemento. 

Esta família define-se em termos de um certo 
número finito de parâmetros, ditos deslocamentos 
generalizados, e é em termos destes desloca- 
mentos generalizados que se escrevem as condi- 
ções reduzidas de compatibilidade que definem o 
modo de interacção dos vários elementos. 

Seja pois q” o vector dos deslocamentos gene- 
ralizados. Os campos de deslocamentos consenti- 
dos no interior do elemento são do tipo 


Ú == ês qó (63) 


em que LÍ é uma matriz cujos elementos depen- 
dem das coordenadas. Para cada vector Qº tem-se 
assim um certo campo de deslocamentos. 

A expressão da energia de deformação de um 
elemento é 


mação: É af 
Uº = 2 /s e dO (64) 
Q 


em que 9º é o subdomínio correspondente ao 
elemento. 
Introduzindo (54) e (55), vem 


Ut= > | (Du)"H (Du) do 


65 
Q (65) 
Introduzindo (63), tem-se finalmente 
e * e | 
U=54 Kg (66) 
em que a matriz simétrica 
e e,”T e 
= Q 

K Ju H(DL)d (67) 


se diz matriz de rigidez do elemento. A matriz 
de rigidez conserva assim a sua função de matriz 
simétrica da forma quadrática da energia de 
deformação. 
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3.4.3 — Exemplo 


Seja por exemplo um elemento triangular com 
campos de deslocamentos lineares, isto é, da 
forma 


u=C,+C, X,+C,X, (68) 
u=C,+C, X, +C; X, 
(68) é equivalente a 
1 E ME DIDO Tas 
: Et XX, | (69) 


em que Cº é o vector dos coeficientes C,, 
Cu + o og 

Tomem-se agora para deslocamentos genera- 
lizados os deslocamentos dos vértices, isto é, 
(ver fig. 13), 


f b 
a 
Fig. 13 
|] 
qt= |? 
a 
Us 
E 
14 (70) 
| 
c 
ui, 
c 
U, | 
Em virtude de (65), vem 
o e e ne 
q M G =NC 
Me? (71) 
id 


em que as matrizes Mº*, Mº” e MºS se obtêm 
particularizando para os nós a, b,ec a matriz 
Mº (cujos elementos dependem das coordenadas). 

Observa-se que Nº é uma matriz quadrada de 
sexta ordem. 
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Tira-se de (71) 


Cce=nº qº (72) 
pelo que 
u= MN qf (73) 
Comparando com (63), vê-se que 
Le= MN (74) 


Tem-se tudo o que é necessário para calcular 
a matriz de rigidez através da fórmula (67). 


3.4.4 — Condições reduzidas de compatibilidade 


As condições reduzidas de compatibilidade 
definem o modo de interacção dos vários ele- 
mentos. 

Exprimem-se como se disse em termos dos 
deslocamentos generalizados. 

Os deslocamentos generalizados são em regra 
as componentes de deslocamento de certos pontos 
da fronteira dos elementos que se dizem nós ou 
pontos nodais. As condições reduzidas de com- 
patibilidade afirmam que os nós coincidentes de 
elementos adjacentes permanecem coincidentes 
após a deformação. 

Seja por exemplo uma estrutura constituída 
por três elementos (fig. 14): 


eco 
Fig. 14 


Seja y!! o vector dos deslocamentos do nó 1 
considerado como pertencente ao elemento I, etc, 
As condições reduzidas de compatibilidade são 


ut! aê ql o. 


u 
ES a TS qua 
| E q DE q (75) 
u!* e gia iai ut 
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Considere-se agora o vector 


(76) 


O vector q diz-se vector dos deslocamentos gene- 
ralizados da estrutura. As equações (75) podem 
tomar a forma. 


q=Al!q (79) 


001: 


* 


As condições de compatibilidade são pois do 
tipo 
(80) 


exactamente como acontece no caso das estruturas 
reticuladas. 


3.4.5 — Análise do sistema 

As equações 

(81) 
definem as forças generalizadas associadas com 


os deslocamentos q”. 
Tem-se pois 
Ut = q" Ksgt = 030º (82) 


Por outro lado, o trabalho das forças exte- 
riores é 


"is 
e à X2ud9+ f poudr (83) 
Q F 
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em que F é a fronteira do domínio 9 correspon- 
dente à estrutura. 
(83) pode tomar a forma 


Ó = FL, Kouao + f p "oudr | (84) 
E NE 
em que Ff F designa a intersecção da fronteira 


de cada elemento com a fronteira da estrutura. 
Introduzindo (63), vem: 


=D] ( inda 14 Led E)3q' = 


E NF 
=3|( fxTL ao + fp” Lear Ata 
q" E nE (85) 
Fazendo 
a LS Da pan | 
F “NE 
(86) 
tem-se pois 
37 =0"5q (87) 
Por outro lado 
js=)U= aU (88) 


Logo 


0T5q=5 0sq'=( 5 0Aº) dq (89) 
e + e f 


Sendo 34 arbitrário, vem 
eT e 
0=5aº q (90) 
e 
que representa as equações de equilíbrio da 
teoria. 
Dispõe-se assim do seguinte conjunto de equa- 


ções inteiramente análogas às que se obtiveram 
para as estruturas reticuladas. 


equilíbrio : 


E 
O= 2º qº (91) 


compatibilidade : q)==A" q (92) 
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forças-deslocamentos : 0º = Kº q“ (93) 


A eliminação dos Q“ (método dos desloca- 
mentos) leva imediatamente, como no caso das 
estruturas reticuladas, a 


Q=K q (94) 


em que 


o à 
K= DA KAS (95) 


3.5 — Aplicação do método aos restantes 
modelos analíticos. 


As equações simbólicas apresentadas em 3.4.1 
são, como se afirmou, válidas em todos os 
modelos analíticos da Teoria das Estruturas. 
Mantêm-se efectivamente na teoria das placas, 
das cascas, das barras, embora mude o signifi- 
cado e o número dos elementos dos vectores 8, 
e, U ep e as expressões dos elementos das 
matrizes H, D, E e N. 

A expressão (67) da matriz de rigidez de um 
elemento e (86) do vector das forças exteriores 
actuantes ma estrutura foram estabelecidas a 
partir das equações simbólicas, pelo que se man- 
têm também em todos os modelos analíticos. 


3.6 — Aplicações às peças lineares 


A matriz de rigidez de uma barra recta de 
secção constante foi determinada em 2.2. Dese- 
ja-se agora determinar a mesma matriz de rigidez 
considerando a barra como elemento finito. 

Hã que discretizar o campo no interior da 
barra, isto é, exprimi-lo em termos dos desloca- 
mentos generalizados, como foi feito em 3.,4,2 
para a elasticidade plana. 

Acontece porém que, enquanto no caso da 
elasticidade plana essa discretização era aproxi- 
mada, no caso da barra o campo de desloca- 
mentos pode exprimir-se rigorosamente em 
termos dos deslocamentos Us e us e das rotações 
9º e 9º das extremidades. 

Tem-se 


Us = & “a x, 6º + & us + %, Be (96) 
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em que 
mid PO a (97) 
E E 
xo x 
x I ; 
q = — 2, + x (98) 
1º L : 
9 
x, 2X, 3x r (99) 
ê E L2 
x x 
=" —— (100) 
L L 


Desde que não se considere o trabalho da barra 
à tracção, o vector dos deslocamentos generaliza- 
dos da barra é 


| ue] 
= | 
“ES | ; (101) 


pelo que, comparando (96) com (63), se conclui 
que 


Le= [2 &o &s 2, | (102) 
A energia de deformação é 
L 
2Ut= [ (uspetax = 
0 
l. E] 
eta?” “1 | (a, &a a, 2, | dx, qó (103) 


em que as linhas denotam derivações em ordem 
AX 

Conclui-se de (103) que a matriz de rigidez 
(matriz da forma quadrática da energia de defor- 
mação) é: 


L L L 
Kº=EI [ (ot)? dx, | af a; dx, [ a) a; dk, [af a; dx, 

= 

Rs aê EA a ço "d 

| aj q; dx, | (07) dx; [07 05 dx; | of a; dx, 
0 | ) “o 

mL L E al 

| as a, dx, | as as dx, | tos)! dx, )s a; dx; 
“0 “0 

a SL sL er à aê 

| aj o, dx, | EM a; dx, Re a dx, J (a) dx, 
“0 "0 Lo 

(104) 
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Introduzindo (97), (98), (99) e (100), obtém-se 


| 


Ke=| 125; ses = ad 
L L L pr | 
| we qa 8 go 
| E" L L L 
mass Es 125 SigEs 
| L L L E? 
| | | 
6E1 EI GEL EL 
L L r E! 
(105) 


São exactamente os mesmos coeficientes de 
rigidez que tinham sido obtidos em 2.2. 

Quanto às forças exteriores, dado que a estru- 
tura se reduz à barra, a sua expressão é 


0 = | Lº xdO + | LC pdF (106) 
o" ENE 
Supondo que a carga é uniformemente distri- 


buída, de densidade p, X = p. 
Tem-se pois 


E Ski 
Re 
poa | p 
à la, 12 E. 
na A Did o! o” 
o x, | P 5 
|—p— 
| aa |) 


As forças generalizadas exteriores são pois 
agora as reacções e momentos de encastramento, 
Lembra-se que se chegou à mesma conclusão em 
2.9, raciocinando com forças de fixação. 

Conclui-se pois que, em virtude de a equação 
(96) ser exacta, as fórmulas (67) e (86) condu- 
zem a expressões exactas para a matriz de rigi- 
dez e forças exteriores. Exaetos são portanto ainda 
os resultados finais da análise fornecidos pelo 
método dos elementos finitos, o qual de resto 
não é mais, repete-se, que uma generalização dos 
métodos clássicos de análise das estruturas reti- 
culadas, 
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Os resultados deixam de ser exactos desde 
que se saia fora do domínio da teoria das peças 
lineares, pois que a equação (96) deixa de ser 
exacta, isto é, a discretização só se consegue por 
intermédio de aproximações. 


3.7 — Critérios gerais 


O tipo de elemento considerado para a elasti- 
cidade plana está longe de ser o único possível. 
Podem: ter-se por exemplo elementos rectangulares 
ou elementos triangulares dentro do qual o campo 
tenha expressão quadrática e não linear. 

O critério mais importante para decidir sobre 
a admissibilidade de um determinado tipo de ele- 
mento é o da convergência para a solução exacta: 
os valores dos deslocamentos e tensões obtidos 
a partir de análises baseadas em elementos de 
dimensões sucessivamente decrescentes devem 
tender para a solução exacta quando essas dimen- 
sões tenderem para zero. 

A convergência para a solução exacta obriga 
a que a equação de discretização (63) tenha uma 
certa forma. Demonstra-se que, essa forma deve 
ser tal que todos os estados de deformação ho- 
mogénea sejam possíveis dentro do elemento, 
assim como todos os movimentos de corpo rígido. 
No caso da elasticidade plana, isso obriga a que 
as expressões das várias componentes de deslo- 
camento contenham um termo constante arbitrá- 
rio e termos lineares afectados de coeficientes 
também arbitrários. É o que acontece no exemplo 
apresentado em 3.4.3, em que as expressões dos 
deslocamentos só contêm esses termos. No caso 
da teoria das placas, em que as deformações são 
as curvaturas (segundas derivadas do desloca- 
mento transversal), a expressão do deslocamento 
transversal deve conter todos esses termos e ainda 
os termos quadráticos também multiplicados por 
coeficientes arbitrários. 

Outro critério que tem sido utilizado é o da 
compatibilidade. Tem-se procurado que as condi- 
ções reduzidas de compatibilidade impliquem a 
compatibilidade completa, isto é, a continuidade 
dos deslocamentos através das fronteiras dos 
elementos. Isso consegue-se no caso do exemplo 
apresentado em 3.4.3, pois que, dependendo os 
deslocamentos linearmente das coordenadas, man- 
têm-se rectos os lados dos triângulos. Basta por- 
tanto garantir a igualdade dos deslocamentos de 
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cada par de nós coincidentes de elementos adja- 
centes (condições reduzidas de compatibilidade) 
para que os lados comuns dos elementos adja- 
centes se mantenham colados após a deformação. 


4 — CONCLUSÕES 


A Teoria de Estruturas moderna apresenta-se 
livre dos defeitos que caractérizavam a teoria 
antiga. 

Para um engenheiro, interessa sobretudo ter 
bem presente que existe um método susceptível 
de tratar qualquer tipo de estrutura, desde que 
disponha de programas de cálculo automático 
adequados. Em presença desses programas, a ta- 
refa do engenheiro reduz-se essencialmente a 
subdividir a estrutura em elementos do tipo cata- 
logado e em interpretar os resultados fornecidos 
pelo computador. 


BIBLIOGRAFIA 


AsrLunD, 5. O. — Structural Mechanics : Classical and 
Matrix Methods — Prentice-Hall, trand, Co,, Inc. 
Juntando o antigo e o recente, esta obra sueca inclui, 
além disso, apreciáveis introduções históricas dos 
vários assuntos tratados. 

Harr, À, S.; Woonuem, R. W — Frame Analysis, 

Uma das melhores obras dedicadas à análise matri- 
cial de estruturas reticuladas. Notável ainda pelos 
magníficos exercícios que acompanham o texto. 


Livescey, Re K. — Matrix Methods of Structural Analy- 
sis — Pergamon, 1964. 
O livro para quem queira iniciar-se na análise ma- 
tricial das estruturas reticuladas. 


Oviverra, E. R. Arantes — Introdução à Teoria das Es- 
truturas — Lisboa, Ed, do autor, 1966. 
Uma sintese da Teoria das Estruturas tornada pos- 
sível pelas recentes contribuições. 


Pezeminteki, J.S. — Theory of Matrix Structural Ana- 
lysis — Mc Graw-Hill 1968. 
Este livro constitui, com o de Zienkiewicz, uma das 
duas únicas obras que apresentam uma exposição 
do método dos elementos finitos. 


Zrenkiwicz, O. C; Creusa, Y. K. — The Finite Element 

Method in Structural and Continuum Mechanics — 
— Mc Graw-Hill, 1967. 
O primeiro livro totalmente dedicado ao método dos 
elementos finitos. Súmula da grande experiência dos 
seus autores, é, apesar de deficiente no plano teó- 
rico, a obra fundamental para quem queira estudar 
elementos finitos. 


TECNICA Nº 395 


C D. DU. 624,345.611 


NOVA CONCEPÇÃO DO ESQUEMA EQUIVALENTE 


DE UM DIELÉCTRICO SÓLIDO 


Por 


RESUMO 


Critica-se o esquema equivalente clássico usado no 
estudo das descargas parciais e elabora-se fisicamente 
um mecanismo de descarga que conduz ao novo esquema 
equivalente, o qual permite calcular com exactidão as 
grandesas caracteristicas do fenómeno das descargas 
parciais, 


1. INTRODUÇÃO 


O envelhecimento dos materiais isolantes por 
processos eléctricos tem como causa fundamental 
a heterogeneidade dos dieléctricos, através de 
reduzidas porções de materiais estranhos no seu 
seio, que provocam descargas eléctricas intermi- 
tentes designadas por descargas parciais. Em 
geral nos isolantes constituídos por camadas, que 
se usam nos equipamentos eléctricos, formam-se 
cavidades (fig. 1) preenchidas por ar onde, 
a partir de determinado valor da tensão aplicada, 
se desenvolvem descargas eléctricas que vão 
deteriorando o respectivo dieléctrico ao longo do 
tempo. 


i , 1 NE z 4 Ê z co o q . 
E a : o x es A : à | A a : 4 E e 
, “ A z ar 2 Ã A 4; É . “e 
Ea a a] », " ' Rd F, F 2 ra E, a 4 X x % a ES 
ho m h 4 3 “ k 5 e, ' Ed dl El ” je “ as. é a, s 
NONE G ESSE NS SN 
” é . ' À i a! 4, - ad o N E é E 
na e 4 h ] A ; a“ q =. ] a a ho k mo ii 
k ' F, A ' j 
| 


Fig. 1 — Dieléctrico com uma cavidade: região 1— 
cavidade, região 2 — dieléctrico em série com a cavi- 
dade; região 4 — dieléctrico em paralelo com 1 e 2. 


Têm sido publicados diversos trabalhos sobre 
a análise matemática do fenómeno quase sempre 
a partir do chamado esquema simplificado, 
embora por vezes se apresentem esquemas 
eléctricos diferentes para esclarecer qualitativa- 
mente específicos detalhes do fenómeno físico. 
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SYNOPSIS 


The classic diagram used in studies of partial dis- 
charges in a void is criticized. À new diagram is propo- 
sed, which is based on a physical mechanism of 
discharges and which makes possible the accurate calcula- 
tion of the main characteristics of corona discharges. 


O novo esquema que se propõe com base na 
interpretação física das descargas parciais conduz 
a expressões teóricas que satisfazem aos resul- 
tados experimentais divulgados. 


2. ESQUEMA SIMPLIFICADO 


Em 1928 W. Burstyn [1] aó estudar as perdas 
dieléctricas em materiais isolantes constituídos 
por camadas esboçou pela primeira vez um 
esquema equivalente para análise das descargas 
parciais, que lhe permitiu justificar qualitativa- 
mente o aumento do factor de dissipação devido 
ao aparecimento desse fenómeno quando o campo 
eléctrico na cavidade excede um certo valor. No 
seu raciocínio esteve presente a aplicabilidade 
das leis de descarga eléctrica no ar entre eléc- 
trodos metálicos, nomeadamente a lei de Paschen, 
as descargas no meio gasoso das cavidades 
inseridas no dieléctrico, conclusão que em 1919 
se revela no trabalho de F. Dubsky [2]. 

O esquema equivalente (fig. 2-b) é formado 
pela capacidade C, da cavidade e pela capacidade 
complementar C, em série relativa à capacidade 
do dieléctrico entre a cavidade e os eléctrodos. 
Em paralelo com a capacidade da cavidade um 
disruptor representa a existência de descargas 
quando a tensão aos seus terminais atinge o valor 
de disrupção. 

Contudo já em 1922 Delafield du Bois [3] se 
serviu nos seus cálculos de um esquema que cons- 
titui o primeiro passo para o esquema anterior- 
mente descrito não atendendo ao disruptor isto 
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é, abstraindo da existência das descargas parciais 
(fig. 2-a). 

A. Gemant e W. v. Philippoff em 1932 [4] 
definiram o chamado esquema simplificado, já 
hoje designado de clássico na análise das descar- 
gas parciais, pela consideração da capacidade 
C, relativa à restante parte do dieléctrico sem cavi- 
dades em paralelo com a série de C, e C, (fig. 2-c). 


Fig. 2 — Esquema simplificado 
a) Segundo Delafield du Bois (1922); 
b) Devido a Burstyn (1928); 
c) Conforme Gemant e v. Philippoff (19932). 


Este esquema tem servido de base para ilus- 
trar o fenómeno e calcular o número de descar- 
gas por segundo (frequência de repetição de des- 
cargas), a intensidade de corrente e a potência 
aparentes (que se manifestam no circuito exterior). 

A investigação das descargas parciais em die- 
léctricos sólidos, em especial com cavidades arti- 
ficiais, tem contudo revelado particularidades que 
este esquema não esclarece, o que lhe confere 
a designação de esquema simplificado. 

Imediatamente se observa que o esquema não 
atende às dissipações devidas à condutividade 
e não considera a existência de uma corrente trans- 
versal na superfície da cavidade. Além disso 
constata-se uma contradição entre a experiência 
e a teoria desenvolvida a partir do esquema no 
que se refere ao número de descargas por período 
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da tensão aplicada : em geral este número é supe- 
rior ao calculado teóricamente a partir do esquema. 
O facto é consequência de uma segunda contra- 
dição: a tensão da cavidade correspondente 
ao apagamento de cada descarga não é forçosa- 
mente nula, conforme o esquema simplificado 
exige. 


3. INFLUÊNCIA DA CONDUTIVIDADE 


No estudo do fenómeno da ionização nos die- 
léctricos conceberam em 1926 C. L. Dawes 
e P. L. Hoover [5] um esquema que inclui a 
influência das resistências (R, e R,) relativas 
à condutividade, para explicar as perdas die- 
léctricas, em derivação com as capacidades da 
cavidade e do dieléctrico (fig. 3-a), sem a con- 
cepção do disruptor em paralelo com C.. 

Na literatura americana encontra-se em 1958 
o raciocínio de E. C. Rogers [6] no cálculo da 
tensão de inserção das descargas (fig. 3-b) onde 
desprezou o efeito da resistência do dieléctrico 
e deu relevo à condutividade da cavidade. 
Observa-se todavia que a descarga curto-circuito 
R, eliminando-lhe o efeito. 

O aparecimento de maior número de descargas 
que o indicado pelo esquema equivalente simpli- 
ficado foi explicado em 1965 por B. Heller 
e J. Chládek [7] através da acção da resistência 
R, do volume do dieléctrico (fig. 3-c) sobre 
o condensador €, quando a tensão aumenta, pro- 
vocando uma subida da tensão em C, e por con- 
sequência um incremento do número de descargas. 

Ensaios revelaram a possibilidade de surgirem 
apenas duas descargas durante um período da 
tensão. Porém o mínimo valor possível indicado 
pelo esquema simplificado (para V-==2V, na 
expressão (1) a seguir apresentada) é de quatro. 
No sentido de dar explicação a esta ocorrência 
idealizaram em 1963 A. Veverka e A. Hon [8] 
um esquema (fig. 3-d) em que intervém o efeito 
da resistência superficial entre as superfícies da 
cavidade (R |) e a resistência interna do dieléctrico 
entre as superfícies da cavidade e os eléctrodos 
(R,), tendo interpretado matemáticamente este 
efeito como uma alteração do andamento da 
tensão aplicada à cavidade durante os intervalos 
de tempo sem descarga, que em vez de sinusoi- 
dal pura terá uma componente exponencial cuja 
constante de tempo é definida por R,, R,,C,e 
C,; assim quando uma descarga se verifica por 
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Fig. 3— Esquema equivalente considerando o efeito da condutividade 


a) No raciocínio de €C. L, Dawes e Hoover (1926); 
by No pensamento de E. C. Rogers (1958); 

c) Devido a B. Heller e J. Chládek (1965); 

d) Por A, Veverka e A. Hon (1963); 

e) Segundo V, Modoni (1962). 


exemplo no instante do pico positivo da tensão 
aplicada, a tensão na cavidade decai instantânea- 
mente do valor de ignição até ao valor de apa- 
gamento (por mecanismo não explicado) e depois 
sofre um deslocamento tal que vai ocupar apro- 
ximadamente o andamento da tensão sinusoidal 
que existiria se não houvesse descarga, verifi- 
cando-se portanto à segunda descarga no pico 
negativo seguida de idêntico processo de evolução. 

Além das resistências de condutividade finita, 
conforme os elementos capacitivos à que estão 
ligados, indicou em 1962 V. Módoni [9] ainda a 
resistência transversal R relativa à corrente 
transversal sobre a superfície da cavidade (fig. 3-e). 

Basicamente porém em todos estes esquemas 
está o esquema clássico. 


4. NÚMERO DE DESCARGAS POR PERÍODO 
De acordo com o esquema simplificado a tensão 
de extinção das descargas na cavidade terá o 


valor zero, correspondente à descarga total de €, 
sobre o plasma formado, pois fisicamente uma 
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vez iniciada a disrupção aos terminais de um 
condensador, se não intervêm outros factores, o 
condensador descarregar-se-á totalmente. 

Em 1963 W. Boening [LO] admitiu que além 
da cavidade também uma certa zona do dieléc- 
trico à sua volta participa no processo de des- 
carga (fig. 4-a). Uma capacidade C, em paralelo 
com C, por intermédio de uma resistência R, 
representa o efeito da vizinhança dieléctrica e 
uma capacidade C. traduz a restante parte do 
dieléctrico em série com C,. A resistência R, 
corresponde à resistência do canal de descarga. 
À corrente transversal que se manifesta durante 
a descarga percorre a. resistência R, em conse- 
quência da diferença de potencial entre C, e Cy 
quer dizer, quando não houver descarga C, e C, 
têm aos terminais a mesma tensão e a corrénte 
através de R, é nula, 

Com este esquema desenvolvido (fig. 4-b) 
Boening foi conduzido a uma tensão de restabe- 
lecimento V, diferente de zero, devido à acção de 
C, sobre C,, para a qual tende a tensão aos ter- 
minais de C, (fig. 4-c) depois de se anular (con- 
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dição de apagamento da descarga), num intervalo 
de tempo insignificante em relação ao intervalo 
sem descarga. 


Fig. 4 — Esquema equivalente considerando a efeito do 
meio ambiente do dieléctrico, na concepção de Boe- 
ning (1963) 

a) Zona de influência do dieléctrico sobre as des- 
cargas ; 

b) Esquema desenvolvido; 

c) Evolução da tensão aos terminais da cavidade 
durante e imediatamente após a descarga; 

d) Número de descargas por período em função da 
tensão aplicada. 


O pensamento de Boening teve origem no 
desacordo experimental do número N de descar- 
gas por periodo da tensão aplicada e o valor 
teórico previsto pelo esquema simplificado, que 
se deduz obedecer à equação 
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N==4, —s=— (1) 


onde V, é a tensão de inserção (valor eficaz ou 
de pico, conforme o valor de V), ou seja, a ten- 
são aplicada para a qual se iniciam as descargas 
na cavidade. 

Após vários ensaios com diversos materiais 
isolantes onde realizou cavidades de diferentes 
dimensões concluiu experimentalmente que a fór- 
mula deverá ser 


4 V-V, 
IN es teses (cota (2) 
1-p VY 


onde p é igual à relação entre as tensões de res- 
tabelecimento e de ignição das descargas 


(3) 


A expressão (2) degenera na expressão (1) 
para p=0, isto é, para V,/==0, quer dizer, 
quando não houver efeito de vizinhança (fig. 4-d). 

Realizando depois montagens segundo o seu 
esquema desenvolvido constatou para a função 
N (V) andamentos em harmonia com (2). Deste 
resultado salientamos de interessante que existe 
uma influência do ambiente circundante da zona 
em descarga. 


5. MECANISMO DAS DESCARGAS 


Hoje ainda não se encontra esclarecida a com- 
plexidade que constitui o processo das descargas 
parciais. Baseado porém nalguns resultados expe- 
rimentais poder-se-ão elaborar as seguintes supo- 
sições : 


1.4 As descargas verificam-se em zonas restritas da 
cavidade: métodos de sensibilização fotográfica 
(fig. 5) têm mostrado que as descargas se desen- 
rolam em pequenas áreas [11], de preferência 
junto das superfícies longitudinais da cavidade, 


o que é explicado pela menor resistência apre- 
sentada (efeito da condutividade superficial) ; 
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Fig. 5 — Descarga numa cavidade de 5 mm de diâmetro 
e o3 mm de profundidade em politene [11] 


2.4 O meio ambiente exerce acção sobre a evolução 
das descargas: os iões positivos originados pelas 
descargas acumulam-se junto das superfícies da 
cavidade donde resulta uma carga espacial que 
irá participar no fenômeno; 

3.2 Não se realizam descargas simultâneas na mesma 
cavidade: quando uma descarga se desenvolve 
constitui um curto-circuito na zona onde aparece, 
pelo que a tensão se reduz em toda a cavidade, 
impedindo assim o desenvolvimento de outras 
descargas nas restantes zonas da cavidade; 

4.4 Descargas sucessivas surgem em zonas distintas 
da cavidade: após o apagamento da descarga 
numa dada zona a tensão na cavidade eleva-se 
devido ao aumento da tensão aplicada (admitida 
sinusoidal), mas as zonas que não foram sede de 
descarga atingem primeiramente a tensão de igni- 
ção em consequência da distribuição de cargas 
a superfície da cavidade no final da descarga 
anterior ; 

5.4 A distribuição de cargas na superfície da cavi- 
dade é alterada (fig. 6-c): devido ao curto-circuito 
provocado pela descarga as cargas eléctricas 
superficiais, originadas por fenómenos de polari- 
zação, acumuladas das ionizações antecedentes 
causadas por choque dos electrões acelerados 
pelo campo eléctrico aplicado com os átomos do 
gas contido na cavidade ou ainda provocadas 
por fotoionização (fig. 6-b), deslocam-se para 
a zona afectada (fig. 6-a); 
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6.2 Durante a descarga surge uma corrente trans- 
versal na cavidade em consequência da diferença 
de potencial entre a zona de descarga e as 
restantes zonas da cavidade e da existência de 
iões positivos distribuídos pela superfície. 


Ed Ed Ed f Pd * Fi Ed Fa f Fê Pci Pa Pai ade Ed Ed 
Ed Fu = E + ea É: dj E r * 
A a A é 
4 dl ” ES Fa Es A 
rs , ad Pd 
| x | | 
o ERES: E ar | 
b | 
pr ni efron — ano — e me ERES as 
a) 


5) 
E E ar 
o ass bela] su ===" 
, ! 
! “0 X) Pd 
? | —— final E 
Ó naipe e ne . Ee E EE 
| e ao 
-Q(X) E 
| 
Ê 
, e 
Ed 


Fig. 6 — Esquematização das cargas eléctricas 
na cavidade 
a) Cargas superficiais, 
b) Canal de descarga e lotolonização ; 
c) AMeração da distribuição superficial de cargas. 


6. NOVO ESQUEMA EQUIVALENTE 


Com base nas considerações físicas anteriores 
concebe-se um esquema equivalente que traduz 
com aproximação satisfatória a complexidade de 
fenómenos desencadeados no processo das des- 
cargas parciais, que designaremos por esquema 
básico (fig. 7-b). 

A capacidade C,, representa a zona em des- 
carga e C, a restante capacidade da cavidade. 
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As capacidades C, e C€, têm o mesmo significado 
do esquema simplificado. O disruptor simboliza 
de igual modo a realização das descargas. Quando 
a tensão aplicada à cavidade atinge o valor de 
ignição desencadeia-se a descarga, a tensão aos 
terminais de C,, reduz-se rapidamente a zero e a 
descarga termina o processo. Entretanto C, e C, 
descarregam-se sobre o plasma através da resis- 
tência R, na cavidade e transversalmente ao 
campo. Desta maneira a tensão aplicada a C, 
decresce mas não chega ao valor zero, pois quando 
cn atinge a descarga total ainda C, apresenta 
uma certa tensão aos terminais, chamada tensão 
restante em virtude da queda de tensão em R, 
(fig. 7-c). 

Note-se que a tensão restante V, corresponde 
à que Boening denominou por tensão de restabe- 
lecimento mas não tem o mesmo carácter, pois 
as concepções de meio ambiente são distintas 
assim como a evolução da tensão. Para o esquema 
básico V, significa o valor da tensão da restante 
parte da cavidade (exclusão da zona em descarga) 
para o qual a tensão do canal de descarga se 


a) b) 


4) e) 


anula e que portanto resta após o apagamento 
da descarga. 

Quando R,==0 o esquema básico degenera no 
esquema simplificado, sendo então nula a tensão 
restante. Interpretamos este caso dizendo que não 
houve acção do meio ambiente sobre o canal de 
descarga o que, segundo o mecanismo elaborado, 
corresponde à não existência de carga espacial 
apreciável. Assume-se nestas condições que os 
iões positivos originados por cada descarga se 
recombinam com electrões, no intervalo de tempo 
entre duas descargas sucessivas. 

Poder-se-ia conceber C, dividida em duas capa- 
cidades respectivamente em série com C,, e C, 
mas complicaria o esquema sem vantagem pois 
C, é muito menor que C,. De facto C,, pode-se 
desprezar no desenvolvimento matemático perante 
C,. Nestas condições C, é praticamente a capa- 
cidade da cavidade e o esquema toma a forma 
do esquema reduzido (fig. 7-d), onde não se 
indica o disruptor para acentuar que se trata de 
um circuito apenas válido durante o período de 
descarga. A resistência R, em rigor inclui a resis- 


6. | EM l 
dE Cc, —— C —- 3 
pad prece camocalica east obs «o 


- 


ig. 7 — Novo esquema equivalente 


a) Zonas da cavidade (De 1) e do dieléctrico (2 e 3); 


bj Esquema básico ; 


c) Evolução da tensão aplicada à cavidade durante 


a descarga ; 
d) Esquema reduzido ; 
e) Esquema completo. 
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tência do canal de descarga, sendo esta variável 
entre valores elevados no início e no final da des- 
carga e com valores muito baixos durante 
a descarga. 

Finalmente, se se atender às condutividades 
no dieléctrico e na cavidade, representadas pelas 
condutâncias G,, G, e G,, obter-se-á O esquema 
completo (fig. 7-e), que apenas tem valor quali- 
tativo pois a consideração das perdas faz-se atra- 
vés do factor de dissipação. 

Em resumo, o esquema reduzido permite o tra- 
tamento analítico do fenómeno das descargas 
parciais com notável simplicidade. 


7. CONCLUSÕES 


Fez-se a análise crítica dos esquemas equiva- 
lentes referidos na literatura sobre descargas par- 
ciais nos dieléctricos sólidos com cavidades. Atra- 
vés de um mecanismo físico do processo de des- 
carga construíu-se um esquema que responde às 
objecções a que o esquema clássico está sujeito. 

A confirmar a validade das hipóteses efectua- 
das poder-se-á deduzir, a partir do esquema pro- 
posto, a fórmula (2) induzida experimentalmente 
para o cálculo do número de descargas por 
período e calcular-se a intensidade de corrente 
e respectiva potência em perfeito acordo com os 
resultados dos ensaios. 
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C. D. U. 624.341.5/94 


ELEMENTOS SOBRE A PRODUÇÃO E O CONSUMO DE ENERGIA 
NA REDE ELÉCTRICA NACIONAL (*) 


— Breve nota mensal 


Quanto às afluências, há a referir que as «afluências 
internas», por efeito das reduzidas precipitações, atin- 
giram apenas 26º/, da respectiva média; o caudal 
transitado em Miranda apresentou-se inferior ao caudal 
natural, por efeito de a albufeira de Ricobayo ter pas- 
sado a armazenar, vindo o regime de afluências ao 
Douro a situar-se em cerca de 66º, da média; as 
afluências ao conjunto do sistema produtor foram ape- 
nas de cerca de 366 GWh, valor que corresponde a 
uma produtibilidade hidroeléctrica de 0,46 e tem uma 
probabilidade de ser excedido da ordem dos 88 “/, 

Em face das afluências desfavoráveis, intensificou-se 
novamente o ritmo de desarmazenamento no conjunto 
das albufeiras; o desarmazenamento durante o mês de 
Dezembro foi de 163 GWh, situando-se o enchimento, 
no final do mês, em cerca de 764 GWh, valor este que 
corresponde a cerca de 36º/, do máximo possível, 


II — Elementos gerais (GWh) 


a) Mensais Variação 

1968 | 1969 | “Vo 
Produção hidráulica Didi na O E 527,9 ai 
Produção térmica (P+)..... “6,1 |— 
Produção total (PT)...... as 538'8 614,0 + 10 


Energia recebida de = 


não pertencentesao RNC (Er) | 53 2585 — 51 


Exportações (Ex)......,..... 1,8 2,2 |-+- 22 
Importações (1)... «ccsece.s. 4,1 24,4 +-495 
Saldo importador (Sp)....... + 2,5 | +22,2 /+865 
Cons. em bomb, hidroelé. Ch) (*)| 18,9) 00 —100 
Produção para con- | | (1) 
sumos perman. (Pcp)..... 498,9 | 585,4 417,8 
Produção “para con- | | 
sumos não perman, (Pcnp) .. 23,6 53,1 — 1 
po a 066,4 | 638,5 + 13 
Coeficiente de hidraulicidade 1,22 0,46, — 


NOTA: 

(1) O aumento percentual da produção para consumos permanen- 
tes, tendo em conta a incidência dos domingos e dias especiais, 
é respectivamente de 19,1 e 10,7º/g. 


HI — Diagramas de carga dos dias característicos 


4. feira: 

1968 o 1969 

Produção hidráulica  (P ) MWh 20028 17 27% 
Produção térmica (Pe) MW | 1400 | 4 620 
Produção total (Pr) MWh [21428 (2189 
Trocas com | Export. (lx) MWb 205 50 
Espanha À Import. (1 )MWh | 20 | 50 
Cons. em bomb, hidroel, (Cb) MWh| 1280 U 
Prod, p.* cons. perm. (Pep) MWh 18 043 21021 
Prod.p.*cons.não perm.(Penp)MWh| 1920 1655 


TOTAL 19963 [226% 


e — mes mo 


P T+(- Ex) MWh. 


Potênciamax. MW 1 poi l Ser 

eu x, Potência min. MW 490 
Ê : Pr Os] Utiliz. da ponta horas 19,0 18,7 
ERA Factor de carga 0,79 0,78 
“a —|— - = e 
Us Potência máx. MW, 968 1142 
Es p Potência min. MW 404 457 
Um cp | Utiliz, da ponta horas 18,6 18,4 

| Factor de carga 0, 7 0, 17 
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DEZEMBRO 


IV — Energia armazenada nas principais albufeiras 


No fim do mês 


Albufeiras: | pm 
GWh 0% 43) 


CD me - -. — ——— | 


Alto Rabagão . ..... «| 462,2 47,5 
Paradela, ». cc vc co os 68,4 30,1 
vanda Novas saw «| 106,9 83,5 
Balamonde . « «vcs co és 25,6 92,8 
Caniçada «eso soe ss 6,0 19,6 
MAMÊ ss DIS ER ka Ea 36,6 36,4 
Cabril ; ssa ria go é a|' VASO 33,9 
Castelo do Bode +... +. 90,1 30,1 
Guilhofrei . .. gue a a 2,5 80,1 
Lagoa Comprida . . . «(2156 40,3 
Santa Lisias é & 5 os é . 27,8 45,1 
Pracana é» : cs o . 0 0 60 5,4 41,9 
POVOS o vam ie wc ué (3) 4,6 32,1 
o" com A, Rabagão . . .| 997,1 43,7 
Potal | sem A, Rabagão. . +. 464,9 40,5 
NOTAS 


(1) Coeficiente de enchimento em energia das albufeiras. 


(2) Inclui 1,5 GWh armazenados em Vale do Rossim no início do 
mês e 0,7 GWh no fim do més, 


(3) Inclui 1,2 GWh armazenados no açude do Poio no início do mês 
e oo GWh no fim do mês. 


(*) Elementos extraídos das estatísticas mensais do Repar- 
tidor Nacional de Cargas (R. N, C.). As produções e 
os consumos das empresas do R.N.C. representam 
cerca de 940/y dos totais do Pais, 
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CONSTRUCOES TECNICAS, LDA. 


PRAÇA DO MUNICÍPIO, 13-3º TELEF. 366506 LISBOA 2 » DELEGAÇÕES: CABO VERDE - ANGOLA - MOÇAMBIQUE 


E 


DIRECÇÃO DE SERVIÇOS DE CAMINHOS DE FERRO, PORTOS E TRANSPORTES (DSCFPT) — 
Moçambique — Porto de Lourenço Marques — Consolidação, Reforço e Aprofundamento do Cais 
Gorjão entre as Carvoeiras. 


TÉCNI CA — XV II O miRa = 25000 
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se! 
et 


A 


ZA 


—" ss 


Si, So qd! uumtms ”, 
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Real” .. 


FEDERAÇÃO NACIONAL DOS PRODUTORES DE TRIGO (FNPT) — Silo de 
Concentração de Lisboa (Beato) — Empreitada global de construção: 3.º fase em 


curso para ampliação da capacidade, de 42 000 t para 60 000t. Fundação com 
estacas QD 1,00 metro. 


CONSTRUCOES TECNICAS, LDA. 


€. D. U. 663.63: 620.193 + 620.191 


O PROBLEMA DA CORROSÃO E DAS INCRUSTAÇÕES NOS 
EVAPORADORES PARA A DESSALGA DA ÁGUA DO MAR 


RESUMO 


Na primeira parte deste artigo analisa-se a experiên- 
cia obtida até agora, no que respeita à corrosão, e pro- 
curam indicar-se as melhores soluções, 

Na segunda parte trata-se o problema das incrusta- 
ções descrevendo os métodos existentes para as evitar e 
finalmente determinam-se as condições para um funcio- 
namento sem incrustações. 


1. A CORROSÃO 
1.1 A experiência obtida 
1.1,1 — Introdução 


Uma fracção importante dos custos da água 
dessalgada é devida ao custo da própria unidade. 
É pois necessário construir esta com materiais 
baratos, mas ao mesmo tempo capazes de durar 
a vida do evaporador, sem necessitar de manu- 
tenção elevada. 

Até agora a maior parte das experiências pro- 
vêm de ensaios laboratoriais, instalações piloto 
e das instalações pará demonstração dos diversos 
processos construídas pelo OSW (Office of 
Saline Water, Department of Interior, U.S.A.) [9]. 

Na última década a maior parte das grandes 


por  J SARAFANA 
Eng.º Químico 1. 8, T. 


SYNOPSIS 


The begenning of this paper aims at providing 
an analysis of corrosion in evaporators and used 
practice. 

The second part looks on the scale formation and 
prevention in the desalination plant and how to operate 
without troubles. 


unidades de dessalga construídas foram do tipo 
MSF (Multi Stage Flash) pelo que os dados 
disponíveis se referem quase, unicamente a estas 
unidades. No entanto, os resultados podem ser 
extrapolados para os outros tipos de evapora- 
dores [11]. 

Vamos pois analisar as diversas partes do 
evaporador MSF para daí concluir dos melhores 
materiais em cada uma delas. 

Para os outros evaporadores estes resultados 
serão uma indicação a ser confirmada experi- 
mentalmente. 


1.1.2 — Tubos de permuta de calor 


De um modo geral pode dizer-se que o funcio- 
namento destes tubos é bom qualquer que seja 
a liga escolhida (quadro 1). 


QUADRO | 


Análise de 55 unidades MSE para uma vida média de 4 anos 


TIPO DE LIGA 


[Aninio! 90/10 | 70/30 | 
Ne de tubos. css ca VR SA O O NR Ss ela É RS E EIN RS | 532.775 | 49.816 | 110.000 7.940 
"o de tubos substituídos +... cce. So e o ii DC É A 0,47 | 0,89 0 o | 
d/, de área substituída +... cc... SE a O RD SST z | 0,38 0,96 0 ) 
Nº de eváporadore + cce das umas Ea ia ER a7 18 16 1 
"de evaporadores reparados +... «cv... a 38 | 44 0 8) 
Vida dos tubos (max. e min.) (em meses) +... «cce rcraas 10-88 | 17-79 24-124 26 
Vida média dos tubos (em meses) ...ciccccactcrccves | 50 | | 
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| Fo | 
| Liga de | Cu-Ni| Cu-Ni| ed 
Titanio | 


48 | 71 26 
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Entre as ligas de cobre, a liga Cu — Ni 70/30 
é a melhor. À liga de alumínio embora pareça 
melhor que a de Cu — Ni 90/10 é preciso ter em 
conta que muitas unidades usam tubos de 
Cu — Ni 90/10 nos andares de maior temperatura. 

Com efeito calculando a percentagem de área 
substituída nas diversas secções teríamos valores 
diferentes (quadro II). 


QUADRO IH 


Vu de área substituída 


" Cu-Ni Liga de 

Tipo de Liga 90-10 alumínio 
Secção de rejeição . . 1,00 “fo 0,52 "q 
Secção de recuperação . .| 0,25 */ 0,34 “Jo 
Secção de aquecimento . 0,93 */y 1,67 “fo 


Esta comparação indica-nos que a liga Cu-Ni 
90/10 é apesar de tudo algo melhor que a liga 
de alumínio nas secções de recuperação de calor 
e de aquecimento |2] [3] [4.. 

A maior velocidade de corrosão nos andares 
de rejeição de calor deve-se a perfurações dos 
tubos causadas por conchas marinhas fixadas nos 
tubos e a resultante corrosão localizada o que 
indica a necessidade de projectar a tomada de 
água do mar de forma a evitar a introdução de 
vida marinha nos tubos. 

Embora o número de substituições seja pequeno, 
há sempre uma corrosão uniforme em todo o 
tubo que depende da liga e das condições de 
funcionamento. 

É de esperar uma velocidade de 0,25-1,52.107 
cm/ano no caso da liga Cu-Ni 20-30, e de 
0,29-1,56.10 * cm/ano no caso da liga de alu- 
mínio. Parece portanto ser de contar com uma 
vida de 20-30 anos de serviço. 


3 


1,1.3 — Armação dos tubos e separadores de vapor 


A armação dos tubos é de uma liga de cobre 
variável enquanto que os separadores de vapor 
são de monal [4]. 

No que diz respeito à deterioração destas peças 
o problema não é grande. O problema de manu- 
tenção mais frequente é o das incrustações que 
se depositam nos separadores. 
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1.1.4 — Canal da água 


A manutenção do canal da água é um dos 
maiores problemas de conservação deste tipo de 
unidades, contribuindo grandemente para os 
períodos de paragem [2] [3] [4] [8]. 

O problema mais vulgar é a erosão e a cor- 
rosão combinadas da superfície do canal, resul- 
tante das mudanças de velocidade e de direcção 
da salmoura quente, 

Tem-se verificado corrosões elevadissimas tais 
como de 1 cm em cada 400 h de funcionamento 
para o ferro fundido ou de 0,3 cm/ano no caso 
de aço. 

A fim de reduzir a corrosão têm sido ensaiadas 
coberturas de plástico (PVC), recobrimentos 
metálicos (à base de cobre) e protecções catódicas 
das ligas de ferro. 

O uso da liga Cu-Ni 90/10 já foi ensaiado 
numa dezena de unidades com relativo êxito. 
O seu preço, não favorece porém a expansão nas 
unidades de maiores dimensões. 

Têm sido feitos estudos no sentido de cons- 
truir este canal em betão armado. Ultimamente 
uma firma americana anunciou ter descoberto 
um processo de misturar o betão com plásticos 
o que conferiria à mistura uma grande resistência 
e impermeabilidade. 


1.1.5 — Às bombas 


Não há ainda unanimidade quanto à melhor 
liga para as bombas usadas no evaporador [4]. 

Usam-se normalmente bombas de bronze e de 
aço inoxidável ou ligado (Ni). 

Parece haver tendência para o uso de bombas 
de aço sempre que haja possibilidade de as 
bombas estarem sujeitas a cavitação. 


1.1.6 — Câmaras de «flash» 


Desde que a salmoura seja desarejada e os 
não condensáveis expulsos o aço é suficiente- 
mente resistente à corrosão para poder ser usado 
nas paredes sem qualquer revestimento, sendo 
suficiente uma revisão anual para realizar uma 
reparação apropriada. No caso de maiores corrosões 
a sua cobertura com plástico não oferece 
problemas [2] [3] [4]. 
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1.1.7 — Razões para quebra de produção 


A capacidade média das unidades consideradas 
foi até agora 11º/o maior que o total das capa- 
cidades de projecto. 

A maior parte das unidades (quase 60/) tinha 
um factor de utilização superior a 90%. As 
razões que obrigaram a parar ou a diminuir a 
produção foram por ordem decrescente de impor- 
tância [4]. 


1) — Limpeza ácida das incrustações 
2) — Manutenção geral 
3) — Manutenção da tomada de água 
4) — Reparação do canal de água 
5) — Substituição dos tubos por entupimento 
6) — Desequilíbrio das temperaturas na unidade 
7) — Desequilíbrio das pressões entre a central 
de produção de energia e o evaporador 
(unidades mistas) 
8) — Problemas com a caldeira (unidades 
simples) 
9) — Ajustamentos dos sistemas de vácuo e 
de purga 
10) — Elevada salinidade do produto 
11) — Avaria das bombas 
12) — Reparação da secção de aquecimento 
13) — Problemas com o funcionamento do 
sistema de tratamento da alimentação 
14) — Falta da procura para a água produzida 
15) — Problemas com o sistema de recirculação 
16) — Problemas nos orifícios de comunicação 
entre os andares e na realização da 
evaporação «flash» por este motivo. 


1.2 — Conclusões 
1.21 — Introdução 


A escolha dos materiais para os evaporadores 
é muito importante e está relacionada com 
factores económicos e técnicos pois pretende-se 
produzir água ao mínimo custo possível. 

É sempre aleatória a previsão das avarias numa 
máquina e o custo da paragem devido a avaria 
ainda é mais difícil de determinar sobretudo 
(caso da água) se se pretende conhecer o prejuizo 
feito à comunidade devido à falta de água. 

Apesar de tudo é sempre possível tecer algumas 
considerações baseadas nos factos experimentais 
que acabámos de ver. 
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1.2.2 — Tubos de permuta de calor 


Se é necessário um serviço absolutamente 
livre de precalços torna-se quase obrigatório o 
uso de tubos de Cu-Ni 70/30 ou de titânio pois 
neste caso o custo mais elevado destes tubos 
pode justificar-se devido a menor manutenção 
necessária e ao diminuto número de paragens 
de produção. 

No caso contrário, podemos usar uma liga 
mais barata (Cu-Ni 90-10 ou uma liga de alumínio) 
o que obriga a maior manutenção e a um mais 
elevado número de paragens mas, diminui imenso 
o preço da primeira instalação de unidade. 

Neste caso o projecto do evaporador deve 
facilitar a substituição dos tubos. 

Atendendo à experiência de que já se dispõe 
é também conveniente o uso da liga de alumínio 
nos andares de rejeição e de recuperação de 
calor e da liga Cu-Ni 90/10 nos andares de 
recuperação de calor e na secção de aquecimento. 

Se os tubos de titânio vierem a baixar de preço 
as grandes produções (construídas para funciona- 
mento contínuo) podem tirar grandes vantagens 
da sua resistência à corrosão pela água do mar, 
pois o facto de não haver práticamente paragens 
para manutenção permite, para a mesma produção 
construir uma só unidade (o abastecimento está 
assegurado com muito maior certeza) e essa 
unidade será mais pequena que uma outra em 
que se tivesse de descontar os tempos de paragem. 


1.2.3 — Armação dos tubos e separadores de vapor 


Uma vez que os materiais já usados não 
apresentam problemas de maior no que diz 
respeito à erosão ou corrosão qualquer deles 
serve. No entanto é de investigar o possivel uso 
de outros materiais mais baratos. 


1.2.4 — Canal de água 


É evidente que o problema da manutenção 
deste canal vai obrigar a estudar a escolha do 
material e a rever o projecto clássico deste canal 
de forma a melhorar a economia global. 

O uso de ferro fundido ou aço inoxidável 
diminui os encargos de primeira instalação mas 
obriga a uma protecção adequada. 

O uso dum novo material, combinando o 
cimento convencional e um plástico, agora em 
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estudo pelo OS W, capaz de resistir à erosão, ao 
choque térmico e à corrosão e além disso imper- 
meável à água e ao ar, parece vir a ser prome- 
tedor. 

Actualmente se se pretende uma instalação 
com um elevado factor de utilização (>> 95"/,) 
deve usar-se uma estrutura barata (ferro fun- 
dido ou aço) recoberta com uma liga à base de 
cobre capaz de resistir à corrosão. 


1.2.5 — Às bombas 


A experiência recolhida até agora permite 
apenas rejeitar as bombas de ferro fundido para 
a circulação de água salgada e a escolha de pás 
de aço inox para maior resistência à erosão por 
cavitação e à corrosão. 


12.6 — Câmaras de flash 


As câmaras dos evaporadores MSF estão 
sujeitos a uma corrosão muito pequena, a não 
ser que o desarejamento seja mal feito ou que 
a temperatura da salmoura seja elevada. O reves- 
timento com um polímero resolve normalmente 
este problema nos casos mais agudos. Também 
se usam revestimentos de materiais como o vidro 
ou o cimento em algumas partes das câmaras de 
vaporização. 


1.2.7 — Controlo da corrosão 


A fim de controlar o estado do evaporador e 
de avaliar as causas da corrosão a fim de poderem 
ser ao menos atenuadas deve instalar-se nos 
evaporadores aparelhagem capaz de medir em 
qualquer instante: 


a)— O grau de polarização da corrosão elec- 
troquimica. 


h) — As perdas de metal através de medidas de 
resistência eléctrica. 


:) — A espessura dos tubos por meio de ultra- 
sons ou da atenuação de intensidade de 
um feixe de raios É ou y. 


d) — À temperatura, o pH e o valor do oxigé- 
nio, anidrido carbónico e ácido sulfídrico 
dissolvidos, em vários pontos do evapo- 
rador. 
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Sempre que a unidade pare e também para as 
peças substituídas devem fazer-se os exames: 


a) — Perda de massa das peças substituídas. 
b) — Exame visual da corrosão. 
c) — Análise dos produtos da corrosão. 


d) — Análise da vida média dos tubos. 


128. —- À alimentação 


A vida dos materiais de construção do evapo- 
rador em especial e dos tubos de permuta de 
calor é muito influenciada pela presença de con- 
chas e de vida marítima na água de alimentação. 

O estudo duma nova unidade deve pois incluir 
um projecto para a tomada de água de um sis- 
tema que permita remover eficientemente tudo o 
que possa vir a alojar-se nos tubos e causar a 
corrosão localizada destes. 

São também de avaliar as variações sazonais 
da composição da água do mar sobretudo naquelas 
substâncias que podem vir a afectar o rendimento 
do evaporador tais como a poluição com sulfu- 
retos, migrações de molúsculos e acumulação de 
algas marinhas. 


2 — AS INCRUSTAÇÕES 
2.1 — Introdução 


Chamam-se incrustações aos depósitos duros e 
aderentes que geralmente se formam nas super- 
fícies de troca de calor e que interferem com o 
bom rendimento da operação contínua. 

Os meios que dispomos para o controlo das 
incrustações dependem das variáveis operatórias 
como temperatura máxima da salmoura, concen- 
tração máxima da salmoura e quantidade de sal- 
moura recirculada. 

Estas variáveis operatórias e a eficiência do 
evaporador estão tão estreitamente ligadas com a 
economia da conversão que justificam todos os 
esforços feitos no sentido de precisar com seégu- 
rança os limites de formação das incrutações. 

A sua formação deve-se à inversão da solubi- 
lidade a quente de certas substâncias. 
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2.2 — Os primeiros passos 


Os primeiros estudos feitos procuraram resol- 
ver o problema das incrustações de carbonato de 
cálcio e de hidróxido de magnésio (Fig. 1) [9] 
[10] [11] [12] [13] [14] [15] [16] [17] [18]. 


alcolindade 
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factor de concentração 


120 
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a SO 
[= 
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pH 
Fig. 1 


Com efeito a água do mar está potencialmente 
saturada em carbonato de cálcio bastando para 
isso decompor térmicamente o ião bicarbonato 
em solução. 


2HCO 7 — CO;*+CO,+H,O 
O cálcio combina-se com o ião carbonato for- 
mado, formando-se precipitado de CaCO, (ara- 
gonite). 
Cat? + CO, * —» Ca CO, 


Acima de 70ºC dá-se a hidrólise parcial do ião 
carbonato: 


Co;”* + H,O — CO, + 2 OH” 
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o que faz aumentar o pH da solução fazendo 
precipitar o magnésio presente (brucite) 
Mg"* + 20H” — Mg (OH), 

Qualquer destes constituintes das incrustações 
alcalinas torna-se menos solúvel à medida que a 
temperatura sobe pelo que as precipitações ten- 
dem a ocorrer e a ser maiores no filme de tran- 
sição entre o líquido e a superfície de transfe- 
rência de calor. No entanto, se aumentarmos a 
pressão do anidrido carbónico estabelece-se o 
equilíbrio a maior temperatura e portanto a 
reacção de decomposição do Ca CO, não se dá 
enquanto essa temperatura não for atingida. 

Também as paredes da superfície de permuta 
de calor podem facilitar ou dificultar a formação 
de incrustações na medida em que a sua rugosi- 
dade facilite a formação dos núcleos de germinação. 
A corrosão destas superfícies vai pois fazer 
aumentar as incrustações. 

Para evitar a formação das incrustações alca- 
linas têm-se ensaiado vários métodos tendo 
começado por ser usada a adição de uma mis- 
tura designada por Hagevap ou P. D. 8 e com- 
posta por trifosfato de sódio, derivados sulfónicos 
de lenhine e um agente anti-espuma. 

Os polifosfatos têm a propriedade de evitar a 
deposição do Ca CO, mesmo em soluções sobre- 
saturadas quando adicionados em quantidades 
muito pequenas da ordem de 1-5 p.p.m. Porém 
acima de 90ºC a degradação térmica dos polifos- 
fatos torna-se excessiva devido à velocidade de 
degradação aumentar (o seu valor dobra cada vez 
que a temperatura aumenta de 6'C) formando-se 
um fosfato menos polimerizado e por último o 
ortofosfato. 

A eficiência com que este tratamento funciona é 
muito variável. 

Com efeito muitas vezes a quantidade de lamas 
formadas é tão grande que obriga a uma lavagem 
com ácido demasiado frequente. 

A falta dum conhecimento fundamental sobre 
o mecanismo do tratamento não permite melho- 
rá-lo com facilidade pois só se pode fazer inves- 
tigação por tentativas. 

Existem, no entanto, estudos sobre a influência 
do ferro, da sílica e do cloro da água de alimen- 
tação, na eficiência do tratamento. 

Também tem sido ensaiado o uso de polímeros 
orgânicos mas por enquanto nenhum é tão bom 
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como polifosfato para o tratamento das incrus- 
tações alcalinas [19] [20] [21] 


2.3 — O controlo das incrustações pelo ajus- 
tamento do pH 


Vimos anteriormente que as incrustações alca- 
linas eram o resultado da decomposição térmica 
dos bicabornatos. 

A adição de ácido forte (normalmente H,5SO, 
por ser o mais barato) vai fazer libertar o ani- 
drido carbónico antes deste passar a carbo- 
nato [9] [10] [11] [12] [13][14] [15] [16] [17] [18] 
[22] [23] 


HCO, 4H” — CO, + H,O 


Normalmente a edição de 100-120 p.p.m. de 
ácido sulfúrico corresponde à neutralização 
de 809 do bicarbonato o que é suficiente para 
evitar as incrustações. Antes do evaporador a 
água de alimentação é desgaseificada. 

Como é natural, é de evitar a adição de ácido 
em excesso pois daí só resulta maior custo ope- 
ratório devido ao gasto de reagente e à maior 
corrosão do evaporador. 

Se o tratamento se fizer em boas condições 
o pH à saída do desgaseificador deve andar 
entre 7,2-7,6. 

Este tratamento é muito mais seguro que o 
Hagevap e custa aproximadamente o mesmo desde 


20 40 BO BO 100 


que o ácido necessário não seja encarecido pelo 
transporte. 

Por outro lado podemos aumentar a tempera- 
tura máxima até ao limite de solubilidade do sol- 
fato de cálcio (122 ºC) o que permite uma maior 
economia térmica do evaporador (Fig. 2). 

Note-se que as linhas a cheio indicadas na 
figura referem-se ao equilíbrio. A cinética do pro- 
cesso de precipitação é porém muito importante. 
Com efeito, podemos ultrapassar a linha de equi- 
líbrio da anidrite sem que haja qualquer precipi- 
tação desde que a taxa de aquecimento ou de 
evaporação sejam baixas. Na figura indicam-se 
para um e outro caso as percentagens de sulfato 
de cálcio precipitado em função da concentração 
relativa e da temperatura [24]. 


2.4 — Os novos processos em desenvolvimento 
2.41 — Electrólise da salmoura 


Devido ao elevado custo do transporte de ácido 
pensou-se em realizar a electrólise da salmoura 
como fonte do ião hidrogénio sempre que as dis- 
ponibilidades de ácido fossem reduzidas [25]. 

A salmoura que sai do evaporador segue em 
parte para a célula electrolítica. No primeiro 
andar desta célula usam-se eléctrodos de aço 
macio para produzir iões de ferro que vão ser 
oxidados no segundo andar da mesma célula 
onde se usam eléctrodos de grafite (fig. 3). 


- um de solub 


por evap (ºlo) 


por aquec (lo) 


temperatura (ºC) 


Fig. 2 
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Fig. 3 


No entanto, o processo é muito caro não só 
devido ao consumo de electricidade como as custo 
e vida dos eléctrodos, à capacidade da célula e à 
vida do separador entre os dois andares. 

Com efeito a baixa utilização dos ânodos e o 
entupimento do separador devido a incrustações 
dos hidróxidos de magnésio e de ferro formados 
por hidrólise encarecem demasiado o processo. 


2.4.2 — Adição de anidrido carbônico 


A adição de anidrido carbónico sob pressão 
foi já proposta como um meio para evitar a for- 
mação de incrustações alcalinas. 

O aumento da pressão CO, evita a reacção de 
decomposição dos bicarbonatos evitando a forma- 
ção de carbonato de cálcio [26]. 

A água de alimentação será desgaseificada 
após a passagem na secção de aquecimento e antes 
de chegar ao evaporador própriamente dito, 
a fim de limitar a corrosão deste pelo oxigênio 
dissolvido. 

O anidrido carbónico é separado e recirculado 
e o carbonato de cálcio aí formado e precipitado 
é removido para evitar o entupimento e a erosão 
do evaporador. 

O processo no entanto não parece prometer 
devido aos maiores custos da bombagem a ele- 
vada pressão e à dificuldade duma separação efi- 
ciente do Ca CO, sem uma perda de carga ele- 
vada. 


2.4,3 — Recirculação de núcleos de cristalização 
Esta técnica consiste na edição de núcleos de 


cristalização à água de alimentação com o intuito 
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de fornecer uma superfície para deposição das 
incrustações [27]. 

Com efeito é de esperar que se se adicionarem 
núcleos de Ca CO, e de Mg (OH), a precipita- 
ção se faça mais facilmente sobre estes que sobre 
as paredes do evaporador. 

A salmoura terá no fim que ser filtrada e os 
finos recirculados para o evaporador. 

A tentativa de usar núcleos de Ca SO, e subir 
a temperatura para além dos 121ºC não resultou 
pois obriga a um controlo muito apertado da 
granulometria, da velocidade de crescimento dos 
núcleos e da estrutura dos núcleos para que estes 
não se depositem, formando uma espécie de 
cimento nas zonas estagnadas (figura 4). 
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Fig. 4 


Além disto, era necessário adicionar o Ca SO, 
nos andares quentes do evaporador pois de outro 
modo ele dissolvia-se para vir a precipitar mais 
tarde, nas paredes de transmissão de calor. 

Apesar das dificuldades o método foi aperfei- 
çoado, [28] e espera-se que a introdução de uma 
coluna de leito fluidizado a seguir à secção de 
recuperação de calor e na qual estaria contida 
a secção de aquecimento virá resolver o problema 
das incrustações de Ca SO, (fig. 4). Actualmente 
o maior problema será o de diminuir a perda de 
carga na coluna de tratamento. 
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2.4.4 — Resinas permutadoras 


A permuta dos iões cálcicos por iões de sódio 
foi estudada com o fim de possibilitar um maior 
rendimento do evaporador [29] [30]. O seu maior 
inconveniente reside porém nos elevados capitais 
de que este processo necessita para poder tratar 
a água de entrada duma unidade de dessalga 
industrial. 

Além disso o custo da resina necessária para 
substituir aquela que se degrada com o tempo é 
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2.4.5 — Qutros procéssos 


Foram ainda estudados processos de remoção 
dos constituintes das incrustações por precipita- 
ção e por osmose inversa ou osmose. 

Este processo requer evidentemente grandes 
áreas de membrana o que indica uma unidade 
industrial para o que o desenvolvimento actual 
ainda não é suficiente [32]. | 

O primeiro (fig. 6) parece atractivo do ponto 
de vista económico mas exige ainda um aperfeiçoa- 


subprodutos 


a ca 


alimentação 


NaOH) Ss ou Ha SO 


áqua tratada 


Fig. 5 


também importante pelo que os custos operató- 
rios não podem ser baixos. 

Saliente-se o estudo técnico-económico duma 
unidade deste tipo feito pela Catalytic Cons- 
truction Co em colaboração com o Office of 
Saline Water e cujo diagrama simplificado se apre- 
senta na figura 5 [31]. 

Este relatório prevê para uma unidade de des- 
salga produzindo do 200.10 m /dia a instalação 
duma unidade com a capacidade de tratamento 
de 400.10” m'/dia. 
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mento mais completo para minimizar o custo dos 
produtos químicos usados [33]. 

Qualquer destes processos envolve capitais 
elevados pelo que a sua economia é favorecida 
no caso de grandes unidades. 


2.5 — Condições para um funcionamento sem 
incrustações 


No actual estado de desenvolvimento do 
controlo das incrustações nos evaporadores, o 
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Este simples ensaio 
mostra como reduzir 
OS CUSTOS 
de manutencão 
dos sistemas hidráulicos. 


Mesmo que a sua fábrica seja do tipo médio, tem 
grandes probabilidades de estar a dispender mais 
do que deveria em reparações dos sistemas 
hidráulicos e também em consumo de lubrificantes. 


Foi por isso que criâmos este ensaio simples. Ele 
não serve apenas para revelar problemas na sua insla- 
lação, mas também para lhe dizer com antecipação 
quanto poderá economizar, resolvendo esses mesmos 
problemas. 


Tudo isto faz parte do PROGRAMA MOPBIL PARA 
REDUÇÃO DE CUSTOS aplicado aos sistemas hidráu- 
licos - estudado especificamente para lhe proporcionar 
economias na manutenção. 


Na realidade, após estudo das suas instalações e 
na sua própria fábrica, estamos aptos a dizer-lhe que 
economia, em escudos, lhe poderemos proporcionar. 
E isto é apenas em lubrificantes e custos de repa- 
Este equipamento exclusivo, de que faz parte uma seringa para ração. Aquilo que economizará na produção, prove- 


verificar contaminações nos úleos, pode ajudar os engenheiros 


da Mobil à requer o consumo de « leo tudraulicos na sua fabrica niente das operações hidráulicas, e um bónus. 


Mobil 


vendemos mais vendendo menos 


TÉCNICA = XVIII] 


